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Abstract 

We propose a discrete fractional random transform based on a generalization of 
the discrete fractional Fourier transform with an intrinsic randomness. Such dis¬ 
crete fractional random transform inheres excellent mathematical properties of the 
fractional Fourier transform along with some fantastic features of its own. As a pri¬ 
mary application, the discrete fractional random transform has been used for image 
encryption and decryption. 
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1 Introduction 


It is well known that the mathematical transforms from time (or space) to 
frequency domain or joint time-frequency domain, such as Fourier transform, 
Winger distribution function, wavelet transform and more recent fractional 
Fourier transform, etc. have long been powerful mathematical tools in physics 
and information processing. For instance, Fourier transform has been the ba¬ 
sic tool for signal representation, analysis and processing, image processing 
and pattern recognition. In physics the Fourier transform describes well the 
Fraunhofer (far held) diffraction of light and thus has been the fundamental of 
information optics [1]. More recently, in the research of quantum information, 
Fourier transform algorithm has been adopted as an effective and fundamental 
algorithm in quantum computer [2]. Wavelet transform is a kind of windowed 
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Fourier transform (Gabor transform), however with variable size of the win¬ 
dows [3]. Therefore wavelet transform has been an extremely powerful tool 
in signal representation in time-frequency joint domain with multi-resolution 
capability, and has been extensively used in image compression, segmentation, 
fusion and optical pattern recognitions. 

The signihcance of mathematical transforms manifests itself further when the 
fractional Fourier transform was re-invented in 1980’s [4,5] and became ac¬ 
tive from 1990’s after its physical interpretations were found in optics [6,7,8]. 
Actually at the beginning, Namias tried to solve Schrodinger equations in 
quantum mechanics using fractional Fourier transform as a tool with little 
notice in community. However, the fractional Fourier transform has found it¬ 
self in real physical processes of light propagation in a graded index (GRIN) 
hber, which is equivalent to a near held diffraction of light, Fresnel diffraction, 
with a quadratic factor. Thus the fractional Fourier transform can be easily 
realized in a bulk optical setup consisting of lenses. The fractional Fourier 
transform provides various of new mathematical operations which are useful 
in the held of optical information processing. And because fractional Fourier 
transform also is a kind of time-frequency joint representation of a signal [9], 
it has found extensive applications in signal and image processing [10]. 

The discrete forms of mathematical transforms have been extremely useful 
in applications, especially in signal processing and image manipulations. In 
fact, discrete transforms can approximate their continuous versions with high 
precision, meanwhile with high computation speed and lower complexities. 
Needless to say. Discrete Fourier transform (or FFT) and discrete wavelet 
transform have been widely used in diherent kinds of applications. Recently, 
discrete fractional Fourier transform (DFrFT) and the relevant discrete frac¬ 
tional cosine transform (DFrGT) have been proposed [11,12]. We have used 
this fast algorithm of fractional Fourier transform in the numerical simulations 
of image encryption and optical security [13,14]. 

As we have demonstrated that the extension of fractional Fourier transform 
have many diherent kinds of dehnitions according to how we fractionalize the 
Fourier transform [15], the DFrFT may also have diherent kinds of versions. 
In our researches of optical image encryption, we ask naturally the question, is 
there any possibility that the DFrFT be random? We have been motivated in 
searching such a random transform because then the image encryption process 
could be simplihed by a single step of transform. Recently we found that, from 
the generalization of DFrFT, we can construct a discrete fractional random 
transform (DFRNT) with an inherent randomness. We demonstrate that such 
DFRNT has excellent mathematical properties as the fractional Fourier trans¬ 
forms have. And moreover it has some fantastic features of its own. We have 
also demonstrated that the DFRNT is a very efficient tool in digital image 
encryption and decryption with a very high speed of computation. The open 
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questions left are concerning with the physical analogies of DFRNT and its 
further applications, which we are considering now. 

We discuss the mathematical definition and properties of DFRNT and provide 
numerical simulation results of the DFRNT’s for one-dimensional and two- 
dimensional signals in the following sections in details. 


2 Mathematics of the discrete fractional random transform 


We begin our discussions from the dehnition of DFrFT proposed by Pei et al 
[11]. A one-dimensional DFrFT can be expressed as a matrix-vector multipli¬ 
cation 

X,(n) = F“x(n), (1) 


where x(n) is the inpnt vector which has N elements, F" is the kernel trans¬ 
form matrix and a is the fractional order. When a = 1, the DFrFT becomes 
the DFT as X( n) = Fx(n), with matrix F indicating the kernel matrix of 
DFT. 

The transform matrix is dehned as follows. Firstly, because the fractional 
Fonrier transform has the same eigenfunctions with the Fonrier transform, we 
can calculate the eigenvectors {Vj} [j = 1,2,, N) of DFrFT with a real 
transform matrix S of discrete Fourier transform which is defined as follows 
[16] 


S 


2 1 0 0 
1 2 cos cn 1 0 

0 1 2 cos 2a; 1 


0 

0 

0 


1 

0 

0 


( 2 ) 


1 0 


0 0 ... 1 2 cos(A^- IV 


where u = 2 'k/N. The matrix S actually is not the kernel transform matrix 
F of the discrete Fourier transform (DFT). However, the matrix S commutes 
with matrix F, i.e. SF = FS. Thus the eigenvectors of S are also the eigenvec¬ 
tors of F, only they correspond to different eigenvalues. Because the matrix S is 
symmetrical, the eigenvectors {Vj} are all real and orthonormal to each other. 
They form an orthonormal basis which is equivalent to the Hermite-Gaussian 
polynomials in the case of continnous Fonrier transform and fractional Fonrier 
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transform. Therefore, from the matrix S we obtain the N x N eigenvector ba¬ 
sis matrix V, which is formed by the N column eigenvectors {vi, V2, ..., v^r} 
as V = [vi V2 ... Vjv]. In the calculation of DFrFT, the DFT-shifted version 
of eigenvectors {Vj} are taken. 

Next step, we determine the eigenvalues of DFrFT. We know that the eigen¬ 
values of the continuous fractional Fourier transform can be written as 

Afc = exp(—iQ;A;7r/2), /c = 0,1, 2,..., cx). (3) 


In DFrFT the eigenvalues are not changed, only we have a limited numbers of 
eigenvalues being taken into account, say k = 0,1,2,N. Those eigenvalues 
again construct a matrix Dq, as follows 


J diag (^1, e ..., e jy jg 

I diag (^1, . . . , e-*a(^-2)7r/2^ ^-iaNnl 2 'j jy jg gyg^^ 


where Dq, is an TV x TV diagonal matrix. It must be noted that in the ex¬ 
pression of Dq, there is a jump in the last eigenvalue for the TV is an even 
integer. Such assignments of eigenvalues of DFrFT are consistent with the 
DFT’s multiplicity rules [17]. 

So now we have already had the eigenvector basis V and the corresponding 
eigenvalues Dq. In the hnal step, we can express a kernel transform of DFrFT 
by the eigen-decomposition. The transform matrix of the DFrFT can then be 
defined as 

F“ = VDqV'^, (5) 


where indicates the transpose of the matrix V. Because the eigenvectors 
are orthonormal, we have VV^ = I and I is the identity matrix. And also 
D_q = D* . Those relations conform that the DFrFT have the same properties 
with the continuous FrFT. 

The above procedure gives a formal way to construct the DFrFT. The matrix 
S is the most important figure in the both DFrFT and DFT, because from 
it the eigenvectors of DFT and DFrFT can be calculated. The matrix S is 
symmetrical and therefore its eigenvalues are all real and the eigenvectors are 
orthonormal. The transform kernel matrix of DFT and DFrFT can thus be 
constructed by eigen-decomposition with different eigenvalues. 

How does a discrete fractional Fourier transform becomes random? The essence 
of the generalization from the DFrFT to DFRNT is to change the matrix S to 
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a random matrix. The DFRNT can be defined by a symmetric random matrix 
Q. The matrix Q is generated by an iV x iV real random matrix P with a 
relation of 

Q = (P + P^)/2, (6) 


where we have Qik = Qki- Similar to the definition of DFrFT, we can also gen¬ 
erate N real orthogonal eigenvectors • • •) of matrix Q. Those 

eigenvectors can be normalized by the Schmidt standard normalization pro¬ 
cedure. Then we have N orthonormal eigenvectors {v/ji, Vki 2 i ■ ■ ■ ■, From 

those column vectors a matrix 

^R = [vri Vr2 . . . Vrrt] (7) 

can be achieved, where rV'ii = I- We do not need to DFT-shift the eigen¬ 
vectors here, because these eigenvectors are generated from a symmetrical 
random matrix. 

The coefficient matrix, that corresponds to the eigenvalues of DFRNT can be 
defined as 


, , 27rQ; \ / Ana \ 

Dr« = diag ( 1,exp ( 1 ,exp J 


,exp 


2{N — l)7rQ;' 
M 


( 8 ) 


which is used in the process of eigen-decomposition. One should notice that 
the eigenvalues of DFRNT are not necessarily relevant to those of DFrFT’s 
or DFT’s. However, we take a similar form so that the DFRNT may have 
similar mathematical properties. In Eq. (8) there is no jump for odd and 
even integer N. We introduce here an integer number M in the coefficients. 
It indicates the periodicity of DFRNT with respect to the fractional order 
a whose significance will be shown bellow. The kernel transform matrix of 
DFRNT can thus be expressed as 

R“ = VrDr^V^. (9) 


Therefore the DFRNT of a one-dimensional discrete signal is written as 


XR,„|(n) = R"x(n) 


( 10 ) 


The expansion of DFRNT for two dimensional signal is straightforward as 

Xr(q,) = R"x (R")^. 
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The most important feature of DFRNT is that its transform kernel is random, 
which results from the randomness of matrix Q, so that the result of transform 
is totally random. The eigenvectors of DFRNT depends upon the random 
matrix Q (or P), therefore if we change the matrix P, then the results of 
DFRNT’s is different. Furthermore, the DFRNT inheres most of the excellent 
properties of DFT (or DFrFT). It can be easily verihed that the DFRNT has 
the following mathematical properties. 


Linearity, DFRNT is a linear transform, he. R“(ax + 6y) = aR“x + 6R“y, 
where a and b are constants. 

Unitarity, DFRNT is a unitary transform, i.e. R“" = (R")* because 
= D|j^. The inverse DFRNT exists and is dehned as R “. 
Additivity, DFRNT obeys the additive role as the DFrFT (and FrFT) 
does for the fact that R“R^ = R^R“ = R"+^. 

Multiplicity, DFRNT has a periodicity of M in our dehnition, i.e. R"+^ = 
R", where we can change this periodicity with changing the integer M. 
Parseval, DFRNT satisfy the Parseval energy conservation theorem, i.e. 
Y^N-1 I -vA / _ sr^N—1 


Xa(R)ik) =Em=o\xi^)\ 


What we are interested most is randomness of DFRNT and the information 
retrieval capability because of its multiplicity. As the fractional order a = 
IM, where I is another integer, the DFRNT output signal ^R{a) return to 
its original function x. Otherwise, it will be totally random when the order 
a 7^ IM even a small aberration occurs. While with the order a = /M/2, i.e. 
at the half of its period, the output signal is real. Such a fascinating 

feature can be rigorously proved in mathematics, because in this case Drq, 
is real (and therefore the kernel transform matrix R" is real). It may also 
be intuitively seen if one recall that the Fourier transform has the property 
of = f{—x). The only difference is that the amplitude of DFRNT, 

when a = /M/2, is random but not |x(—n)|. 


To illustrate the basic feature of DFRNT, and to make a comparison with 
DFrFT, we present here the simulation results for one-dimensional rectangu¬ 
lar signal using DFrFT and DFRNT, in Fig. 1 to Fig. 3, respectively. The 
signal has a rectangular window with period of [40,60]. The total number of 
points is A = 100. The numerical results of DFrFT is given in Fig. 1, with 
the fractional order a = 0.25, 0.50, 0.75 and 1.00, respectively. We can see 
that the amplitudes of DFrFT’s gradually change from the original signal to 
its Fourier transform. However, the results will be different for the cases of 
DFRNT. We calculate DFRNT with two different formats of random num¬ 
bers, one is normally distributed random number (illustrated in Fig. 2) and 
the other is uniform random number within [0,1] (Fig. 3), as the matrix P. 
Here we set the periodicity of DFRNT as M = 1. From the results shown in 
Fig. 2 and Fig. 3, one can see that the amplitudes are all randomly distributed 
when the fractional order a ^ 1. Whereas when a = 1, the output signal goes 
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back to its original, i.e. the signal recovers. 


Note that when a = 0.50, the phase values are 0, tt, or —tt whatever random 
matrix is used. That means we always get real transform results when the 
fractional order a is the half of the periodicity M. The amplitudes for a = 0.25 
and a = 0.75 are the same, however their phases are conjugated. That is 
because ^R{M-a) = ^R{-a) = . When a = 1, the amplitudes of 

transformed signal totally retrieved, however, there exit about ±7r/2 phase 
fluctuations at the components of x{k) = 0. 


3 Image encryption and decryption: a primary application of DFRNT 


The primary and perhaps the most important application of DFRNT is cryp¬ 
tography and information security. Because the DFRNT itself is random, the 
DFRNT of a two-dimensional signal can be directly used for image encryp¬ 
tion with any fractional order a ^ IM. The decryption process is simply an 
inverse DFRNT. The main encryption key is the matrix Q. We have known 
that changing a random matrix Q indicates changing the present DFRNT to 
another DFRNT. The results are unrelated even they have the same fractional 
order a. The DFRNT is quite sensitive to the random matrix and also to the 
fractional order, which have been demonstrated in our numerical simulations. 

Numerical results of DFRNT’s of an image are given in Fig. 4 and Fig. 5 with 
normally and uniformly distributed random matrix respectively. The test im¬ 
age is the photo of henna with the size of 256 x 256 pixels, shown in Fig. 4(a) 
and Fig. 5(a). The encryption process is simple, just an a order DFRNT with 
a 7 ^ IM. The encryption results are shown in Fig. 4(b), Fig. 4(c), Fig. 5(b) 
and Fig. 5(c), respectively, with a = 0.50 and a = 0.80. In our simulations we 
let M = 1. The decryption is an inverse DFRNT with same random matrix Q 
and the fractional order —a. However, if we perform the inverse DFRNT with 
another matrix Q which is different from the matrix we use in the encryption 
process, we fail to retrieve the image. The decryption results are shown in 
Fig. 4(d) and Fig. 5(d). And also it is worth to mention that the fractional 
order a can also be served as an extra decryption key. The inverse DFNRT’s 
with wrong fractional orders can not recover the image as shown in Fig. 4(e) 
and Fig. 5(e), with a = —0.502 and a = —0.503, respectively. As we have an¬ 
alyzed from the computation of mean square error (MSE), that the smallest 
security discrimination is about \Aa\min ~ 0.02 for normally distributed ran¬ 
dom matrix and \Aa\min ~ 0.03 for the case of uniformly distributed random 
matrix. When |Aq;| > \Aa\min, one can not visually recognize the image from 
the noisy background. The correct decryption results are shown in Fig. 4(f) 
and Fig. 5(f). The MSE for correct decryption process is in the order of 10“^^ 
for both cases. 
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The image encryption with DFRNT using the fractional order a = /M/2 is of 
great interest for real applications, because then the encrypted image has real 
values. It may be very useful for storing and transferring the secret data in a 
more convenient way, for example, by a photo plate. 

The security strength of the DFRNT encryption can be estimated as 
because the random matrix Q has N{N + l)/2 independent elements. As a 
matter of fact, if one try to search such a random matrix coded with uniformly 
distributed random numbers, the number of steps will be much larger than 
2 N(n+i)/ 2 ^ Therefore, such an image encryption method is considerably secure 
in theory. This encryption algorithm can be easily implemented digitally. 

Is there any applications of DFRNT in Physics? This question now is left for 
an open problem. Nevertheless, DFRNT has a powerful multiplicities with 
changing the formats of matrix Q. The different matrix Q may results in 
different transform with different mathematical properties. Such a feature may 
be helpful in hlter designs in signal and image processing. We can also establish 
an iterative mechanism that the matrix Q can be modihed accordingly. This 
mechanism may be found applications in solving inverse problems in optics, 
such as the phase retrieval. 


4 Conclusions 


We have proposed a new kind of discrete transform, a discrete fractional ran¬ 
dom transform, based on a generalization of the discrete fractional Fourier 
transform. The intrinsic randomness of the discrete fractional random trans¬ 
form origins from a symmetrical random matrix, from which the eigenvectors 
of the transform are generated. The discrete fractional random transform in¬ 
heres the excellent mathematical properties as the fractional Fourier trans¬ 
forms have. A new image encryption and decryption scheme is proposed that 
uses the discrete fractional random transform. Thus the processes of image 
encryption and decryption are very simple, only consisting of an operation of 
discrete fractional random transform and its inverse. Mathematical proper¬ 
ties of the new discrete transform have been given in details with numerical 
demonstration. 
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List of figure captions 


Figure 1. DFrFT of a one-dimensional rectangular window signal. The frac¬ 
tional orders are a = 0.25, a = 0.50, a = 0.75 and a = 1.00, respectively. 


Figure 2. DFRNT of a one-dimensional rectangular window signal with nor¬ 
mally distributed random numbers. The fractional orders are a = 0.25, a = 
0.50, a = 0.75 and a = 1.00, respectively. 


Figure 3. DFRNT of a one-dimensional rectangular window signal with uni¬ 
formly distributed random numbers. The fractional orders are a = 0.25, 
a = 0.50, a = 0.75 and a = 1.00, respectively. 


Figure 4. Numerical results of DFRNT with two-dimensional data using a 
normally distributed random matrix. DFRNT serves as an image encryption 
and decryption algorithm here, (a) The original image, (b) encrypted image 
with a = 0.5, (c) encrypted image with a = 0.8, (d) decryption result for 
image (b) with a different random matrix, (e) decryption result for image (b) 
with a = —0.502, and (f) the correct decryption of the image. 


Figure 5. Numerical results of DFRNT with two-dimensional data using a 
uniformly distributed random matrix, (a) The original image, (b) encrypted 
image with a = 0.5, (c) encrypted image with a = 0.8, (d) decryption result 
for image (b) with a different random matrix, (e) decryption result for image 
(b) with a = —0.503, and (f) the correct decryption of the image. 
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Fig. 1. Numerical simulations of DFrFT of a one-dimensional rectangular window 
signal. The fractional orders are a = 0.25, a = 0.50, a = 0.75 and a = 1.00, 
respectively. 


11 


Phase Phase 




































a=0.25 


a=0.50 


0 20 40 60 80 100 0 20 40 60 80 100 





(d) 


Fig. 2. Numerical simulations of DFRNT of a one-dimensional rectangular win¬ 
dow signal with normally distributed random numbers. The fractional orders are 
a = 0.25, a = 0.50, a = 0.75 and a = 1.00, respectively. 
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Fig. 3. Numerical simulations of DFRNT of a one-dimensional rectangular win¬ 
dow signal with uniformly distributed random numbers. The fractional orders are 
a = 0.25, a = 0.50, a = 0.75 and a = 1.00, respectively. 
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Fig. 4. Numerical results of DFRNT with two-dimensional data using a normally 
distributed random matrix. DFRNT serves as an image encryption and decryption 
algorithm here, (a) The original image, (b) encrypted image with a = 0.5, (c) 
encrypted image with a = 0.8, (d) decryption result for image (b) with a different 
random matrix, (e) decryption result for image (b) with a = —0.502, and (f) the 
correct decryption of the image. 
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Fig. 5. Numerical results of DFRNT with two-dimensional data using a uniformly 
distribnted random matrix. DFRNT serves as an image encryption and decryption 
algorithm here, (a) The original image, (b) encrypted image with a = 0.5, (c) 
encrypted image with a = 0.8, (d) decryption resnlt for image (b) with a different 
random matrix, (e) decryption result for image (b) with a = —0.503, and (f) the 
correct decryption of the image. 
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